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Abstract 


In this article, the Riccati Equation is considered. Various techniques of find- 
ing analytical solutions are explored. Those techniques consist mainly of 
making a change of variable or the use of Differential Transform. It is shown 
that the nonconstant rational functions whose numerator and denominator 
are of degree 1, cannot be solutions to the Riccati equation. Two applications 
of the Riccati equation are discussed. The first one deals with Quantum Me- 
chanics and the second one deal with Physics. 
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1. Introduction 


The Riccati Equation, named after the Italian mathematician Jacopo Fransesco 
Riccati [1], is perhaps one of the simplest and the most interesting first-order 
nonlinear ordinary differential equations. The Riccati equation has the general 


form 
(1) 


where a, band care real value functions, yis the unknown variable and x is the 


y' =c(x)+b(x)y+a(x)y’, 


free variable. 

This equation arises in many fields of mathematics and physics. It has appli- 
cations in random processes, optimal control, diffusion problems, stochastic rea- 
lization theory, robust stabilization, network synthesis [2] [3], and quantum me- 
chanics [4] [5], to name a few. Besides important engineering applications, the 
newer applications include areas such as financial mathematics. 

One of the strengths of the Riccati equation is that it is a unifying link between 


linear quantum mechanics and other fields of physics, such as thermodynamics, 
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and cosmology. For instance, information about the dynamics of a quantum 
mechanical wave packet can be obtained from a complex (nonlinear) NL Riccati 
equation [6]. It has also been shown that complex Riccati equations appear in 
time-dependent quantum mechanics, and the time-independent Schrodinger 
equation can be rewritten as a complex Riccati equation [6]. Moreover, the dy- 
namics of the Bose-Einstein condensate and Friedman-Lemaitre equations can 
be described by a complex Riccati equation [7]. Also by rewriting the solution 
to the equation of motions, which has the Riccati form, and introducing a 
substitution, we can derive well-known expressions from statistical thermody- 
namics [8]. Other interesting characteristics of Riccati’s equation include its 
possession of a superposition formula, the Painlev’e property, and its ability to 
be linearized [9]. 

The Riccati equation also plays an important role in financial mathematics [10] 
since most interest-rate models contain time-dependent functions. In the case of 
Cox-Ingersoll-Ross Interest-Rate Model (CIR), one must find time-dependent pa- 
rameter that ensures equality between model prices and market prices, which in 
finance is known as calibration. Riccati equation is directly relevant to this cali- 
bration procedure, and once closed-formed solutions to a Riccati equation that 
corresponds to the CIR Model are found, these solutions can be utilized to show 
that there exists an analytical expression for any interest-rate derivative in the 
CIR interest-rate mode. For a full-scale reading of the fundamental theories of 
the Riccati equation, see the book of Reid [11]. 

It is well known that there is no general way to analytically solve a Riccati eq- 
uation, and only special cases can be treated, especially if a particular solution is 
known. The solutions of this equation, for instance, can be obtained numerically 
by the forward Euler method and Runga-Kutta. An unconditional stable scheme 
was presented by the author in [12]. An analytic solution of the non-linear Ric- 
cati equation was obtained by El-Tawil et af using Adomian decomposition 
method. Abbasbandy solved one example of the quadratic Riccati differential 
equation (with constant coefficient) by He’s variational iteration method by us- 
ing Adomain’s polynomials [13]. In [14], the authors implemented the Homo- 
topy Analysis Method (HAM) to solve a Ricatti equation. 

The differential transform method (DTM) is a numerical method for solving 
both linear and nonlinear differential equations. It was first introduced by Zhou 
[15]. The method constructs a semi-analytical numerical technique that uses 
Taylor series method for the solution of differential equations in the form of a 
polynomial. Unlike the traditional Taylor series method that takes a long time 
for computation of higher order derivatives, the differential transform method is 
capable of greatly reducing the size of computational work while still accurately 
providing the series solution with fast convergence. 

Moreover, this method can be applied directly to linear and nonlinear ODEs 
without requiring linearization, discretization, or perturbation. 

The article is divided into six chapters. 

In chapter 2, several existing methods of solving the Riccati Equation are ex- 
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plored. In theorem 1, a new integrability condition is given. 
In chapter 3, it is shown that nonconstant rational functions of the form: 


kt+k 
(t) =———._ cannot be solution to the Riccati equation. 
kt +k, 


In chapter 4, the Differential Transform Method is considered. This method 
seems to be very effective and promising for solving Riccati Differential equa- 
tions. 

In chapter 5, applications to Quantum Mechanics and Classical Mechanics are 
discussed. 


Chapter 6 is dedicated to the conclusion. 


2. Solving the Riccati Equation 


It is well known that solutions to the general Riccati equation are not available 
and only special cases can be treated. However, it is always possible to write the 
general solution of the Riccati Equation if a particular solution is known. After a 
change of variable using the particular solution, we obtain the Bernoulli Equa- 


tion, which can be solved by traditional methods. 


2.1. Deriving the General Solution of the Riccati Equation Given a 
Particular Solution 


Consider the general form of the Riccati equation (E): 


&  a(t)y? +B(t)y+C(0) 


Let y, bea solution of (E). That is 


& = a(e)y? + B(t)y, +C(t) 


It is well known that if w=y-—y, then w satisfies the Bernoulli Equation 
(Fi): 
w'+(-2A(t) y, - B(t)) w= A(t)w’. 


which we know how to solve see [16] 


2.2. Methods a Particular Solution to the Riccati Equation Solution 


There are many ways of finding solutions. The techniques are based on a change 


of variable. 


2.2.1. Method 1 
Let (E): y'=A(t)y’+B(t)y+C(t) be the Riccati Equation. We use the 
transformation: 


y=uv-B/A 


where uw is the unknown variable and vis a function of t. We want to find condi- 
tionson A(t), B(t),and C(t) so that (E) is separable equation. 
Substituting yand y’into (E) gives: 
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A’u'v=M —A’u (v'+ Bv) + Au’v? 
where 
M = B'A- BA'+ A’C 


If M=0 and v’'+Bv=0 then the equation is separable, we obtain 


After integrating: 


ejute 
| Avd(c) 


u(t) 


Integrating v’+Bv=0, we have v= e!24) which gives the particular so- 


lution: 
1 


fae Tac) 


u(t) 


2.2.2. Method 2 


We use the transformation 


A C' 
; = 
u’ =(CA)2| 1+ 2A eo 
(cap 
If 
Be AC 
— 5 = constant = K 
(CA)? 


1 
then, u’=(CA)?2 (1 + Ku+ u’) is separable, and we can find an explicit solution. 


2.2.3. Method 3 
Consider the Riccati equation (E): 
“ = A(t)y? +B(t)y+C(t) 
t 

We consider the substitution y(t) =uv-—w, were v and w are the known va- 
riables. 

Taking the derivative of y(t): 

Replacing in (E): 

We obtain 
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uv= -(v' +2A(t)vw+ Bv)u + A(t)w* - Bw+w'+C+ A(t)u°v? 
Assuming that v’+2Avw+ Bv = —bAv*? and Aw*—Bw+w'+C =aAv’, then 
we obtain the following separable equation that we can solve: 
uv= (bu +a+u’) Av’ 


2.2.4. Method 4: Changing the Riccati Equation into a Second Order 
Linear Differential Equations 


Consider the Riccati equation (E): 
y'=A(t)y’ +B(t)y+C 


We can turn the Riccati equation into a second order linear differential equa- 
3 


: . : u 
tion using the transformation y=— " 
u 


, Where wis the new variable. This leads 


to the second order equation: 


ut+{-B- u's acy =0 
A 


C 
We also can use the transformation y= a which leads to the second order 
u 
linear differential equation 
w Cc’ t 
u+ Br u + ACu =0 


' I 


Cc 
If AC is a constant and Bare or eae is a constant, then explicit solu- 


tions to the Riccati Equation can be derived. 


2.2.5. Method 5 
Lemma 1. The Riccati equation (£): 


o = a(t)y’ +b(t)+c(t), where a(t) #0 


can be reduced into the following equation: 


Y _ y? +m(t)y-+n(c) 
dt 
Proot. Using the change of variable y= — and replacing in (E) 


a(t) 
—va(t 


dy _ va(t) va ) 
dt a(t) 


=v" +bv+c(t) > 


a 
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Lemma 2. The Riccati equation (£) can be turned into the equation 
~ =y(t)y’+f(t),where y(t)= Ja(t)}e PO" , where 


t 


wee ai) +¢C 


Y 
Proot. From lemma 1 (E) can be transformed into the equation 


a’ 
var 4[b+Z}yre 
a 
a’ 
Let v--[o+e 
a 


cati equation: 


Theorem 1. /f c(t)= {v0 z e) , then the function 
a 


t)= b(t) a'(t) la(c)| eo Jb(s)at 
y(t) ce a’ ( | a(t) fla(eye mae 


t 
is a solution is of the Riccati equation (E): 


i e()=-[ 50419] (r(ay's where 7(0)=o(0e 


, then the 
function 


is a solution to the Riccati equation (E): 


© = a(t)y? +b(t)y+e(t) 
Proof. According to Lemma (1) the Riccati Equation, 
“ =a(t)y’+b(t)y+c(t) can be turned into (E:) v’=v” (ve Z)vee us- 
a 


ing the change y= z. According to Lemma (2), (Ei) can be turned into (E) 
a 
DOI: 10.4236/am.2022.139049 779 


Applied Mathematics 


M. Ndiaye 


Jeow, where y(t) = ja(t)Je PO. Then w satisfies the Ric- 


M. Ndiaye 


vi=ywet using the change v= -[o 2) + yw , where 
a 


Gr 
a 

Y 
y= ja(t)Je PO". 


, 
t 


Now if c= -[b+2) then wW=yw > =y > wee", Going back- 
a w 


ward we find 


b a’ |a| —[b(t)de 
y(t) (2 + =) a fale "ae 


Same thing if c--[b+=) +y°.Wehave w=yw+y => 
a 


Ma 
w +1 
Jyat= frat => 


arctan w = [rae = 


w= tan({ ydt) 


Going backward we find that: 


y(t)= (2 + “| + el tana tan ({7(t)ae) 


a a 
Example 1 
Consider the Riccati equation 
d 1 1 
i a ar i 
where 
1 


b a! |a| eo Jb(t)de 
aq a 


a flalelae 


and substituting in a(t),b(t) and c(t) gives the solution: 


or ae 
te' e e(t-1) 
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ty t4 
5 2 
A ae t(t—1) 
? : k,t+k, : : 
3. About the Rational Solution y(t)=—*——*. of the Riccati 
k,t+k, 


Equation 


Recall that if the coefficients of the Riccati equation are polynomials, then par- 
ticular polynomials solutions can sometimes be found. For instance, in [17] the 
authors have developed a method by which polynomial solutions to the more 
general Riccati equation: 


A(t)y' = P(t)y? + Q(t)y+R(t) 


can be found, where A(t),P(t),Q(t) and R(t) are polynomial functions in 


t. 
Now how about rational solution in the form 
kt+k 
VO 
kt +k, 


Theorem 2. Let (£): = a(t)y* +b(t)y+c(t) be the Riccati equation. The 


equation (£) is invariant under the rational transformation 


_hh (t)w+y,(t) 
V3 (t)w+y7, (t) 
that is wis a solution of the Riccati equation: 


OW _ 1 (t)w? +1, (t)w+h(t) 
dt 
Proot 

dy (yiw+nw't75)(%3+n)-(mWwtr)(Yw+ yaw + yw) 


2 
2 
= of 2222) +0 HeLa 
Y3wt+ V4 Y3wt V4 


VY sW + NY aW+ VV WH V30 4 FV a NYaW — VW — VoVW—VoVaW' — Vos 


(y,wt+7,) = 
_ hw tmnt (nwtnrwers) awe) 
(y,w+74) (y,w+74) (%3+%) 
(11% — NY) W + (M14 — VF Va%a — Va) Wt (Na —Y2%a) WF V4 — Pos 
=> 


=a(yiw' +2y,7,w+73)+b(n74W + 1%aW+ Yo%aW+ Vora) 
+¢(y,w? +2y,y,w+ 74) 
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(Vrs — Ns) W? + (V4 — Na + 1203 — Va) WH (Na —a¥a) W' +(12%4- 127) 


=(ay? +byyy,+cv3 )w? +(2ay7, + by,%4 + DY 27s + C774) W = 
+(ay; +by7, +72) 
ay, +bY,7¥3+CY3 — N33 
wa 1/3 3 1/3 V3 V2 
NV 4 V2, 
4 LOY 2 + ONY 4 + BY os + 2CP3%4 Ns + Naas t VV, 
VV 4 —V2V3 
4 M2 + DY eV 4+ OV4 = FY 4 t oY 
NV 4 —V213 
; : kt+k, , 
Theorem 3. If the rational function y(t)= Krag 84 20n-constant solu- 
gl + 4 


tion of the Riccati equation © —a(t)y? +b(c)y+e(¢) where k,,k,,k, and 


k, are constant, then the coeftcient a(t), b(t), and c(t) are constant func- 


tions. 
Proot If we choose w=t and k,,k,,k, and k, are constant, then 
kt+k, , . ; ‘ i tate 8 
y(t) oo in the previous theorem. Now replacing w by f¢ in the Riccati 
gl + K, 


equation in the previous theorem, we obtain: 


‘2 a(t)k; +b(t)k,k, +c(t)k; 2 


kk, —k,k, 
4 2a lt) hike +b()(kiks + kak) + 2c(t) kak , 
Kk, ~ kk, 


r a(t)k} +b(t)k,k, +c(t)k; 
k ky —k ky 


for all ¢ 
a(t)k; +b(t)k,k, +c(t)k; =0 
2a(t)k,k, +b(t)(k,ky + kjk, )+2c(t)k,k, =0 
a(t)k? +b(t)k,k, +¢(t)k? = kk, —kyk, 


The Matrix of the system is given by: 


k? k, k, ke, 
M =| 2kk, kk,+k,k, 2k,k, 
k kk, k 


The determinant of Mis given by: 


\t| = kk, +k,k, 2k,k, 7 2k,k,  2k3k, 5|2kk, kk, +k,k, = 
: kk, ki _— k3 ki 5 ks kk, 
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Kk? (Kk ki + kkk; — 2k, kk; ) — kjk, (2k,k,k? — 2k> kk, ) 
+k; (2k k>k, —kykok, —k3k; ) 

= kk? + k?k kjk? —2k?k,kgk? — 2k?k kk? + 2kk2k?k, 
+ 2k, k2k?k, —k,k?k2k, —k3k? 

= kik? —3k?k kjk? +3k,k2k2k, —k23k3 = (kjk, —kyk) 
=|M|=(kk,-kk,) 


Since y(t) is a non-constant solution, then k,k,—k,k, #0. Therefore the 
system has a unique solution (A(t),B(t),C(t)) expressed only in terms of 
k,,k,,k, and k,. Since k,,k,,k,, and k, are constant, then A(t),B(t) and 
C(t) are constant functions. 


4. Differential Transform 


The Differential Transform Method (DTM) is a numerical and analytical me- 
thod for solving a wide variety of differential equations, such as, but not limited 
to, ordinary and partial differential equations. This method provides the solution 
of the given differential equation in terms of convergent series with easily com- 
putable components. The concept of DTM was first introduced by Zhou [4] and 
its main application are in linear and nonlinear initial value problems in elec- 
trical circuit analysis. The DTM gives exact values of the mth derivative of an 
analytic function at a point in terms of known and unknown boundary condi- 
tions in a fast manner. This method constructs, for differential equations, an 
analytical solution in the form of a polynomial. It is different from the Taylor se- 
ries method, which requires symbolic computation of derivatives and involves 
more computations especially for higher order differential equations. Different 
applications of DTM can be found in [18]-[26], to name a few. 

In this chapter, the Differential Transform is used to find explicit solution of 
some Riccati equations. The solutions are expressed in terms of a Taylor series 


expansion. 


4.1. The Essentials of Differential Transform Method 


For convenience, we present a review of the DTM. The differential transform of 
the Ath derivative of a function f(x) is defined as 
a‘ 


Pog eel] 


where f(x) is the original function and F(x) is the transformed function. 


The differential inverse transform of f (x) is defined as, 


(x)= DF (k)(x-%)! 


=0 


From the differential transform and the differential inverse transform we get, 
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showing that the concept of differential transform derives from the Taylor series 
expansion. This method, however, does not evaluate the derivatives symbolical- 
ly. Instead, relative derivatives are calculated by an iterative way which are de- 
scribed by the transformed equation of the original function. 

The fundamental operations performed by differential transform can be rea- 
dily obtained and are listed below in section 4.2. The main steps of DTM is first 
to apply the differential transform to the given equation, which will result in a 
recurrence relation. Then solve the relation and use the differential inverse 


transform to obtain the solution to the problem. 


4.2. Properties 


The following properties can be deduced from the differential transformation 
and the differential inverse transformation and are given below: 

1) If f(x)=g(x)£h(x), then F(k)=G(k)+H(k) 

2) If f(x)=Ag(x), then F(k)=AG(k), where 2 isa constant 

aiff (x)= t then F(k)=(k+1)G(k+1) 


a 
xX 


4) If f(x)= 9) then F (k) =(k+1)(k+2)---(k+m)G(k+m) 
xX 
if k= 
5) If f (x) =x", then F (k) =6(k—m) = : ron 


6) If f(x)=g(x)h(x), then F(k)= 0° G(r) H(k-r) 
7) If f (x)= f(x) f(x): f, (x) then 


F (k)= 3 “SF, (ky) Ey (ky “hE (KK) 
Bx B‘ 
8) If f(x)=e™, then CS aoa) 
9) If f(x)=[og(t)de then (y=), k>1, F(0)=0 


10) If f(x)=sin(@x+a@),then F(k)= o-sin( +c] 


k 
11) If f(x)=cos(@x+a),then F(k) =< cos{ Sa 


k! 
Example 2 
We first start by considering the following equation (E): 
“ =-y' +1, 
with initial condition 
y(0)=0 
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Let Ybe the differential transform of y. 
By taking the differential transform of (E) using properties 1 - 5, the following 


recurrence relation is obtained: 


(K-+1)¥ (K+1) =-Dy (r)¥ (kr) +4(k) 


: i 
Since y(0)=0 then Y(0)= Tey, =0 
If k=0 then Y(1)=->)’ .Y(0)Y(0-r)+6(0) 
Y(1)=-Y(0)Y(0)+1 
Y(ij=1 
Following the same recursive procedure to find Y(2),Y(3),Y(4),¥(5) 
when k=1, k=2,and k=3 wehave: 
If k=1 then 2¥(2)=->.,Y¥(r)¥(1-r)+6(1) 
Y(2)=0 
We find in the same way Y(2) => y(4)=0, Y(5) ==. 
Since x, =0, the differential inverse becomes: 


y= SY (K)x! 


then 
1 2 17 62 
y(t)=t-=t? += t’ + t° 
3 15 315 2835 


Example 3 
Solve (E) 


d 
SY =-e'y? +2ey +e! -e® 


with initial condition 
y=! 
Taking the differential in both sides gives: 


(k+1)¥(k +1) = F,(k,)¥(k, ~k,)¥(k-k,) 
+23 Fy (r)¥ (kr) 4 


0 


(K4a)¥(k+1)=-¥ Fy (4, BY (Ks) 


a 1 (-3) 
29 --y(k 
Fe a 
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The initial condition at x, =0 is: 


1 
v(0)=2 y(0)=1 
If k=0 then 
0 ky 1 0o°o—r 
Y(Q)=-Y SAV (kk )¥(He) +29 ¥ (7) 1+ 
ky =0 k,=0 r=0 
¥(ij=2 
If k=1 then 
4 LZ 1 -3 
2¥(2)=-> VY (k, -k, ¥ (kk, )+2,—¥ (1-1)+ 
kao kao K,! 0 Wr! 1! 1! 
7 
Y(2)=— 
(2)=2 


If k=2 then using the same procedure, find Y(3)= 3 
Substituting Y(0),Y(1),Y(2),Y(3),--- into the differential inverse trans- 
formation, we get the following solution: 


y(t)=14 2x42 -axt 


5. Applications 


In this section, we will see that many equations that arise in physics, cosmology 


and quantum mechanics can also be changed into a Riccati equation. 


5.1. Quantum Mechanics and Schrodinger Equation 


In Classical Mechanics, we describe the motion of an object using Newton’s 
second law which states that the sum of the force (F) acting on an object of Mass 


(m) equals the mass of the object times the acceleration (a) of the object: 
F=ma 


In quantum mechanics that represents the analogy of classical mechanics, but 
at the smallest scale. For example, for electrons orbiting the nucleus of an atom 
Newton’s Mechanics does not apply. Also, in classical mechanics, we describe a 
state of a physical system using the position and the momentum of the objects 
that compose the system. However, in Quantum mechanics, a particle is de- 
scribed using a wave function y(x,y,z), where the square of the magnitude 
lv’ represents the probability that the particle is located at the position 
(x, y,Z, t) , for a single particle moving around a box in three dimension. 

y (% y, z,t) is a solution of the Schrodinger equation: 

ih ow _ h (Gy Oy dy 
2n at 8n°m ax’ ay? a2” 


Jevw 


where /is a constant called then Planck’s constant 
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h 
h = 6.626068x10™ m?-kg/s and h= aA is the reduced Planck’s constant. 
T 


iis an imaginary number i=~-1 


mis the mass of the particle. Vis the potential energy of the particle. 


5.1.1. Relationship between Riccati Equation and Schrodinger Equation 
Let’s consider the one-dimensional time-independent Schrodinger equation 
(one particle) 

aes 
~ 2m @x? 


in Sy (st) -( Jule 


If V(x)=0, we deal with a free motion of the particle. If V(x) = Siw ; 


then we deal with the Harmonic oscillator with constant frequency w=w, or 
time-independent frequency w(t). 

In either cases, the solution of the one dimensional time independent Schro- 
dinger equation: 


fa) hn oO 
ih ,t)=|— 
ae ) 2m oy 


or 


) oo . thors 
ih x,t)=| - +—wrx 
ae Qmox? 2h)” 


can be written as: 
y (x,t) _ A(t) ela(xyt) a A(t Oreo") 


The real-value functions of time a(t), A(t), y(t) satisfy the following 
system of ordinary differential equation 


ts oe (2) 
dt 2h m 
Bh og =0 (3) 
dt m 
dy h 
pir pecs =0 4 
dt om (4) 


(2) is the Riccati equation that we can solve using differential transform. 


Example 4 Choose “ =1, w=e 
Solve: 

—+e"+a’° =0 (5) 
Taking the differential transform in both sides of the equation 


(k1)a(k +t) +Za(ra(k-r) =0, Sa(0)=0 
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If k=0 then a(1)+>.,a(0-r)+1=0, a(1)=-1 
2 =i) a(2)=5=1 using the 


if K=1 then 2a(2)+ Pia (r)a(1-r)+—= 


5 
=-1, 4 Ss 5 Ss = 
a)=—2 as)=——_ a8) — 


same process, @(3) 
The solution of (5) is given by: 


a(t)=-t+t*-t 


5 
oe" t+—t 
6 


5.1.2. Time-Independent Schrodinger Equation and Riccati Equation 


i a 
+(V—-E)w=0 (Harmonic Oscillator 

om pe” tV-E)w=0 ( ) 

Suppose that the particle of mass mis confined in a box that lies in the inter- 


val 0<x<L Let’s consider the case where the potential energy is 


2 
V(x)= TG. 
2 
Let q=mW, so that V(x)= ax The Schrodinger equation then becomes 
in dy (: 2 
+|—x°-E |y =0 
2m dx? 2 r 


Making the change of variable y= ¥  Jeads to the equation 
y 


+ 
4 h? \2 
, 2 My 2m 
yty ee ie 
Solving 
2m 


(k-+1)¥ (K41)+ DY (r) y(k=r) 


r=0 
: —2m —m? 9 
Then we obtain: Y(0)=0, Y(1)= 7 By ¥(ZS0 Y(3)= F*, 
Y(4)=0, 
—4 m 2 3m* i 
5)=——\~E’, Y(6)=0, Y(7)=—— 
E 2 
Choosing w,=— so q=m 
& Wo ; q Re 
We then have: 
—2m m° 23 4m _3 5 3 m* 407 4m 5.9 
= Ex —-— E°x® -——~ E°x® -—— ———E°x’ - 
a nh‘ 5 n° 5H 9 n' 
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y= e 
y 
y' —m—= pn 2S 24 4M 13 6 
dx = | ydx = Ex E°x 
J J a ‘ An‘ 30 n° 
3m* 48 4m 5 10 
40n* 90 n° 


We show that 


2mE mE 
x)= Acos x |cosh] —— x? 


We can find A by Normalizing y(x). We obtain 
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5.2. Applications in Physics—Motion of a Particle in Central Force 


Field 


In a central force field, the force acting on the particle of mass m has the proper- 


ties that: 


The force is always directed from a fixed point towards, or away from the par- 


ticle. 


The magnitude of the force only depends on the distance rfrom 0. 


F is called a central force and the particle is said to move in a central force 


field. The point 0 is reference as the center of the force. 


We have the following properties for a particle moving under the influence of 


a central force. 
1) The path of the particle lies on a plane (planar) 


2) The angular momentum of the particle is conserved 


3) The time rate of change in an area swept by the position vector (from the 


point 0) is constant. 


Using Figure 1, e, =cos@i+sin@j and e, =—sin@i+cos@j, and e, and 


e, are unit vectors. That is, 


from a fixed point 0 towards or away from the particle, then: 
F = f(r)e, 


If f(r) <0, then the force is attracted towards 0. 
If f(r)>0,then F is repulsive away from 0. 


le,/=1 and |le,|=1.Since F is always directed 


5.2.1. Equation of Motion for a Particle under the Influence of a Central 


Force 
The position of the particle is: r=re, 


The velocity of the particle is: v=re, + r6e, 
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Particle 


Figure 1. Illustration. 


The acceleration is: a= (F —r@? Je, + (276 + r6)e, 
Using the Newton’s second law: F = ma 
We obtain 


This is a nonlinear equation that we can turn to a Riccati equation using the 


f 
, r dr 
change of variable u=— where r’=— 


r dg" 


If f (r) = = we, therefore, deal with the Riccati equation 
r 


> 


u'—u? -1=— 
m 


that can be easily solved. 


6. Conclusion 


In this article, a number of integrability conditions of the Riccati equation have 
been discussed in Chapter 2. Method 5 has been added to the list of methods 
that eventually lead to an analytical solution. As the previous methods, the main 
idea is to turn the Riccati equation into a separable Differential equation. For 
each of these methods, the complexity of the method depends on the coefficients 
of the Riccati equation (A(t), B(t),C(t)) : 

The notion of differential transform is also presented in this article. The small 
size of computation and the rapid convergence show that this method is very 
different and could be an alternative for solving non-linear differential equations 


over existing methods. The usefulness of the Riccati equation has been pointed 
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out in the study of the Schrodinger equation in quantum mechanics. Some 
change of variable in the Schrodinger equation leads to the Riccati equation and 
the differential transform is used to write analytic solutions to the Schrodinger 
equation. The motion of a particle under the influence of a central force was ex- 
plored, and again the equation of the motion found was turned into a Riccati 
equation using a change of variable. 

A new result has also been presented in chapter 3 about solutions of the Ric- 


cati equations in the form: 


at+b 
t)= 
( ) ct+d 


where a, b, cand dare constants. It is shown that solutions of this form to the 
Riccati equation cannot exist unless it is a constant solution. A natural extension 


of this article would be to explore the general rational functions of the form): 


where P(t) and Q(t) are polynomial of arbitrary degree. 
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